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ABSTRACT 

The Intent of this project was to develop a course 
for mathematics graduate students at Iowa State University. They 
would design sM write computer programs for use by undergraduate 
mathematics students, and then offer the course and actually produce 
the software. Phase plane graphics for ordinary differential 
equations was selected as the t<^lr-. Prior to the course, the faculty 
coordinators designed the modular structure of the program euid wrote 
some Input/output routines. The course was held, but was plagued by a 
shortage of students. This caused delays, as the program was not 
completed during the course. The software was finally completed 
through a variety of methods, including adapting existing numerical 
program:^, using gradui^te students to write parts of the program as 
master's degree projects, and hiring graduate students to write parts 

the program. The computer program. Phase Plane Graphics for 
Ordinary Differential Equations, is listed in an appendix; it will be 
commerlcally distributed. (Author/MNS) 
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student Produced Advanced Mathematical Software 

Summary 

The intent of this project uias to develop a course for 
mathematics graduate students to design and write computer 
programs for use of undergraduate mathematics studentSi and to 
offer the course and actually ptoduce software The particular 
software project selected was a graphics program for use in an 
undergrc Juate differential equations course. The course was 
helOi but suffered from a shortage of students. Students 
were recruited later to complete the software^ Phase Plane 
Graphics for Ordinary Differential Equations, which will be 
commerc ially distributed. 

Project Director: L.eslie Hogbeni Department of Ma th pinat i c s> 

Iowa State Universityi AmeSi la 50011 
Tel 515-294-8168 or 515-294-1752 

The following materials developed by thi^; project may be obtained 
by calling or writing the atidress above Telephone requests are 
recommendeo for computer diskettes. 

PHASE PLANE GRAPHICS FOR ORDINARY DIFFERENTIAL EQUATIONS 

Input/Output and Graphics Subprograms 

Materials Developed for a Graduate Course In Mathematics 
Software Deve 1 opment 



Project Title Student Produced Advanced Mathematical Software 
Grantee Organization. Iowa State Universityi Ames, la 50011 
Project Director Leslie Hogbeni Department of Mathematicsi 

louia State University* Ames, la 50011 
Tel. 515-294-8168 or 51 5-294--1752 



Executive Summary 



Project Dvervieui 

The intent of this project was to develop a course for 
mathematics graduate students to design and write computer 
programs for use of undergraduate mathematics students, and to 
offer the course and actually produce software. The particular 
software project selected was phase plane graphics for ordinary 
differential equations^ suitable for use in an unde graduate 
differential equations course. Prior to the course* the faculty 
coordinators* Leslie Hogben* Richard Tondra* and Roger Alexander* 
designed the modular structure of the program and wrote some 
input/output routines. The course was held* but was plagued by a 
shortage of students. This caused delays in the project* as the 
program was not completed during the course. The software was 
finally completed through a variety of methods* including 
adapting existing numerical programs* using graduate students to 
write parts of the program as Master's degree projects* and 
hiring graduate students to write parts of the program The 
computer program* Phase Plane Graphics for Ordinary Differential 
Equations* will bi commercially distributed 

Fur p ose 

There is a serious lack of mathematical software at the 
advanced undergraduate ( postcalculus ) level* despite the many 
useful applications of computers to mathematics. There is also a 
lack of integracion of computer use into graduate mathematics 
courses. 

This project attempted to solve both these problems by 
developing a graduate course in which the students would write 
software for use in advanced undergraduate mathematics courses 

Background and Grig ins 

Iowa State University is a large (25*000 students) public 
university offering both undergraduate and graduate degrees The 
Department of Mathematics has more than toO permanent faculty 
members and approximately 90 graduate students Bachelor of 
Science* Master of S lence* and Doctor of Philosophy degrees are 
offered by the department The Iowa State Computation Center 
operates mainframe* minicomputer* and microcomputers* and 
provides consulting services to the university. 

Project Description 

The fall semester of 1984 was spent designing the graduate 
course. Math 5l7x* and making preparations for the software. We 
researched software design and review techniques and decided on a 
format for the cou'^se. Ue developed standard input/output 
subprograms for use in the software. We selected phase plane 
graphics for ordinary differential equations as the program we 



would develop in Spring 1985, ana devised a plan of modular units 
for the program. 

Given a differential equation of the form 

dx/dt = F(x.y) dy/dt O(xiy) 

and an initial points the solution is a path in the Xiy-piane 
There are many applications of such differential equations. 
Because of their importance^ equations of this type are 
frequently studied qualitatively in undergraduate differential 
equations courses^ although quantitative solution of such 
equations is beyond the scope of an undergraduate course. 

The program was to contain the following modules: main 
program (MAIN); to control the region of the x^y-plane that is 
graphed (G'^APH), to enter differential equations (PARSE); to 
graph a single trajectory, i. e , solution, (TRAJECTORY); to graph 
a series of trajectories, i.e. phase plane (PHASE), to graph a 
direction field of vectors (ARROWS); to graph null clines 
(NULLCLN), to find and classify critical points (CRITPT). There 
was also a module to handle menu input/output (MENUS), and for 
technical reasons two small machine-specific modules pertaining 
to graphics and menus 

During the spring semester of 1985, Hogben and Tondra taught 
Math 517x, a graduate course covering numerical methods for 
studying phase planes and developing phase plane software. 
Serious difficulties were encountered with both the quantity and 
quality of students. Only four students registered for the 
course and two subsequently dropped it During Math j17x, two 
riodules were complett?d, MAIN by James Coyle, who subsequently 
drooped the course, and PARSE by Hogben. One student was 
assigned to adapt public domain code in Shampine and Gordan, 
Computer Solution of Ordinary Differential Equations, but did an 
inadequate job. Another was assigned the CRITPT module but did 
not do it 

During the summer of 1985, we hired Coyle to design and 
write the ARROWS module, which he did skillfully and efficiently. 
At the end of Fall semester 1985. Linda Ten Hoeve was given the 
CPITPT module as a creative project required for her Master's 
degree. She completed this project successfully in the summer of 
1986, but did not interface it to the main program. During the 
spring of 1987, Hogben adapted the Ten Hoeve code and Alexander 
adapted the Shamp ine/Gordan code. Practical considerations 
ir^volvmg the use of the program ied to the elimination of the 
PHASE module Thus the program Phase Plane Graphics for Ordinary 
Differential ^quation^> was completed in May 1987 <Z-100 version) 

Because of the difficulties encountered during Math517x, the 
'students did not write the user manual as they went along, as 
originally intended Alexander is currently writing a detailed 
user manual for the program When this is completed, we intend 
to test the software m an undergraduate differential equations 
course In the meantime, brief instructions on how to use the 
progrim are available 

In order to facilitate dissemination of the software, we 
intend to produce versions for other MS-DOS machines, in 
particular the IBM-PC. Technical difficulties have caused 
delaySi but the IBM--PC version should be completed during 1987 
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Project Results 

The principal product of this project is the computer 
program, PHASE PLANE GRAPHICS FOR ORDINARY DIFFERENTIAL EQUATIONS. 
This program is currently available for the ZenitN Z-100 and is 
being adapted for the IBM-PC and similar MS-DOS machines In 
order to obtain widespread distribution* we -are planning to 
distribute it commercially. It has been submitted to CONDUIT 
This program is included in Appendix B 

In the process of deve loping this program* we dev eloped 
menu format input/output and graphics control procedures callable 
from MS-DOS Pascal or FORTRAN. These procedures may be useful to 
other software developers They are included* along with a small 
demonstration program* in Appendix B 

In the process of planning the course* a variety of materials 
wert devel 0 ped These include program design notes* class pi an* 
and a brief software development bibliography These materials 
are included in Appendix C, 

Copies of the ab ove materials may be obtained from 
Leslie Hogben/ Department of Mathematics 
Iowa State University* Ames* lA 50011 
Tel 515-294-8168 or 515-294-1752 

Because information about equipment is necessary to provide 
computer diskettes* such requests are more easily handled by 
telephone. 

Summary and Conclusions 

This project successfully completed the development of 
planned software* Phase Plane Graphics for Ordinary Differential 
Equations. However* the intended method of development by a 
'^.Jass of graduate students was not completely successful. In 
overcoming the student shortage difficulty* two promising methods 
for developing software were utilized Modifying currently 
available research programs* and hiring graduate students. 

Well designed and thouroughly tested programs utilizing 
liopb i st i cated numerical methods are widely available and are 
frequently in the public dofnain. Although such programs 
frequently employ a high level of numerical sophistication* they 
r-ange from difficult to impossible for an inexperienced person to 
use. However* if combined with a user-friendly program to handle 
i npu t/ou tpu t and a we 1 3 wr i t ten manua 1 / existing numerical 
p^^ograms can provide the basis for valuable software. 

The other strategy employed to successfully complete the 
software involved using graduate students to design and write 
parts of the program* but not as part of the class. One student 
wav> hired to write a module and another wrote a module as a 
cTc^ative degree project Both these approaches provide software 
development experience* albeit without the team approach ard 
design review experience originally planned 

The results of this project suggest possible methods of 
developing undergraduate software by incorporating existing 
research subroutines and using gradUdle students, either by 
hiring Ihem or by awarding academic credit for individual 
projects This project also developed a sophisticated and easy 
to use program for use in an undergraduate differential equations 
course* Phase PLane Graphics for Ordinary Di f f er efi t lal Equations. 
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student Produced Advanced Mathematical Software 

Final Report 

Project Overview 

The intent of this project was to develop a course for 
mathematics graduate students to design and write computer 
programs for use of undergraduate mathematics students^ and to 
offer the course and actually produce softaiare The particular 
software project selected was phase piane graphics for ordinary 
differential equationsi suitable for use in an undergraduate 
differential equations course. Prior to the start of the course* 
the faculty coordinators; Leslie Hogben^ Richard Tondra* and 
Roger Alexander/ designed the modular structure or the program 
and wrote input/output routines The course was held* but 
was plagued by a shortage of students This caused delays in the 
project/ as the program was not completed during the course The 
software was finally completed through a variety of methods, 
including adapting existing numerical programs; using graduate 
students to write parts of the program as Master's degree 
pT*ojects> and hiring graduate students to write parts of the 
program. The computi^r program* Phase Plane Graphics for Ordinary 
Differential Equation^?* will be commercially distributed 

Put pose 

There is a serious lack of math emat i cal software at the 
advanced undergraduate ( p os t ca 1 c u 1 u s ) ieveli despite the many 
useful applications of computers to mathematics There 3S also a 
lack of integration of computer use into graduate mathematics 
courses. 
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A substantial body of software exists for eJemnntarg 
mathematics courses (calculus and precalculus) Theie are also 
a large number of sophisticated research level programs 
available^ but thei, are unsuited for use in undergraduate? 
courses Very little software relevant to advanced undergraduate 
(postcalcu lus ) courses is currently available This is primarily 
a consequence of market conditions^ rather than of the 
appropriateness of software to advanced courses While there is 
perhaps a greater ranga o*^ interesting applications in advanced 
mathematics; major textbook publishers currently do not believe 
the development of advanced software is economically viable. 
Once the software is developed* there are means available to 
distribute it. 

The use of programs written by others can make many aspects 
of advanced matheniatics more concrete by allowing the user to see 
a variety of of examples^ much as spreadsheet software can enable 
a businessman to obtain a better understanding of the effects of 
different decisions However^ to thoroughly understand the 
mathematics one needs to understand the und (?r 1 y i n g computational 
procedures (numeric<il algorithms). This is learned by doing the 
programming oneself/ but is beyond the scope of most 
undergraduate mathematics courses It is entirely appropT-iate at 
the graduate level. 

There is a need to incorporate computer programming of 
numerical algorithms into graduate mathematics courses Few 
graduate mathematics courses use computers^ despite the fact that 
niciny graduate students will need to use computers m future 
emp loyment . 
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This project attempted to solve both these problemi:> by 
developing a graduate course in which the students iL:Q;jld write 
softuiare for use in advanced undergraduate mathematics course<i 

Background and Origins 

louia State University is a large (25>000 students) public 
univer«iity offering both undergraduate and graduate degrees The 
Department of Mathematics has more than 60 permanent faculty 
members and approximately 90 graduate students Bachelor of 
Science^ Master of Science^ and Doctor of Philosophy degrees are 
offered by the department The Iowa State Computation Center 
provides academic and research computer services to the 
university. It offers mainframe* miriicomputers and 
mic r ocumputer S/ and provides consulting services. 

The Computation Center was helpful and cooperative 
throughout the project They provided consulting services on 
iTiachine specifics and provided us with some public domain 
graphics software they had written. We had no problems getting 
ici pss to the university's Zenith Z-100 micro computers 

Project Description 

This project was carried out by three faculty members in the 
Department of Mathematics at Iowa State Univorsity> Leslie 
flugben> Roger Alexander* and Richard Tondra* with the assistance 
of several graduate students in the Mathematics Department 

The fall semester of 1984 was spent designing the graduate 
course^ Math 517x; and ma king preparations for the soft ware We 
;e-3earched software design and review techniques and decided on a 
format for the course. We developed standard input/output 
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'^nbprograms for use xn the coftuiare We selected pha«ie plane 
graphics for ordinary differential equations tho pruyram u)e 
liiuld uevelop in GpririQ 19^5, and devised a pl^n of modular un^ts 
for the program 

it was decided that the course »iinuld begin u«ith lectures on 

underlying mathematics that the programmers (graduate 
'itudents) would need for the project. Then each student would be 
assigned a piece of the program to desinn The whole class would 
dji";aiyse the work in design reviews, A design is written in a 
raxture of English and b 1 oc k -s tr u c t u .*ed programming language 
<'5uch as Pascal) Ideally^ a design should be readily 
intelligible to a nonpr ogrammer who is familiar with the 
math emat i c Si and yet should be mechanically translatable into 
programmed code In practice> a design usually starts in 
English and and progresses closer to proniramming language as the 
leview process proceeds The author of the design distributes it 
to the review group (class) and the review group meets at a later 
time to discuss the design The reviewers ask questions to 
clarify the design/ catch cr-^ors/ ask how steps will be 
translated into programming language, etc Materials prepared 
for classroom use in Math 517x, including a sample design, are 
include in Appendix C. 

Phase plane graphics for ordinary differential equations 
was chosen as the software development pioject Given a 
differential equation of the form 

dx/dt - F(x;y) dy/dt = G(x>y) 

^-iod an initial point> the solution is a path m the x^y-plane. 
There are many applications of such differential equations. Such 
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equations are studied in various ways A trajectory (solution 
path) might be needed Critical points (both derivatives 2ero> 
ava studied and classified Because of their importance, 
equations of this type are frequently studied qualitatively in 
undergraduate differential equations coursesi although 
quantitative solution of such equations is beyond the scope of an 
undergraduate coui ie 

The availability of high speed computers has made numerical 
fTjethods (carried out by computer) the best choice for dealing 
uiith such equations. These methods are well developed and 
csiccessible to graduate students Recent developments in 
microcomputers have provided the necessary speed and memory 
capabilities. User convenience combined with superior graphics 
made microcomputers the machines of choice. We had access to the 
University's Z-100 microcomputers. 

During the fall of 1984, the modular structure of the phase 
plane graphics program was designed/ and modules to handle the 
input/output via menus and the graphics window were written 
These modules may be useful to others who wish to design gr^aphics 
^>of|-ware and are included in Appendix B 

The primary considerations in designing the program were 
erise of usei flexibility, and accuracy Accuracy was a major 
lOncern m the implementation of th^ numerical methods In order 
to provide flexibility, the user should be able to enter any data 
lised by the program For example, for graphing a trajectory, the 
initial point, time interval, and error tolerances for the 
numerical routines may be supplied by the user In addition to 
controlling the data for the various graphing actions, the user 
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should be able to change the region of the xiy-plane that is 
graphed and the differential equations 

Ease of use uias a major consideration in the design cf thp 
program as a ui^ole* and the design of input/output. Because 
menus provide t easiest fcrm of user controli it was decided 
that the program should oe menu driven. That iss a list of 
possible choices appears the screen The user selects an 
action from the menu of choices by mov^ing the cursor to the 
desired choice and pressing the RETURN key In the phase plane 
program* the graphics display occupies most of the screen and the 
menu occupies the lower quarter of the screen. To avoid having a 
confusing menu with too many items> there are a number of submenus 
uihich are reached from the main menu Each subsidiary menu 
contains all the options and input data necessary for a single 
action* such as drawing a trajectory. 

Various other features were incorporated into the design to 
facilitate ease of use. Defaults were to be supplied for oil 
data, 50 that a novice user could have the program graph things 
without having to enter data Online helpfiles were to be 
provided. The program was to be crashproof 

The modular structure of the program was designed to 
parallel the menu structure. The program was to contain the 
following modules: main program (MAIN), control of the region of 
the x#y-plane that is graphed (GRAPH), entery of differential 
equations (PARSE); gr^v/oing a single trajectory, i e. * solution, 
(TRAJECTORY), graphing a series of trajectories, i.e. phase plane 
(PHASE)> graphing a direction field of vectors (ARROWS>; graphing 
null clines (NULLCLN), finding and classifying critical points 
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(CRITPT). There (jas also a module to handle menu input/output 
(riENUS)# and for technical reasons tuio 5mall machine-specific 
modules pertaining to graphics and menus. 

On^ other design consideration ujas machine portability. The 
machines uje had available were Zenith Z'-lOO microcomputers. 
Although loosely described as iBM-compa t id 1 e# the graphics and 
screen control on the Z~100 are totally different from the IBM-PC 
In order to make the program available on the IBM-PC and similar 
machines as well as the Z-100# we decided to isolate the machine- 
specific part of the program in small modules (separate from the 
major graphics and menus modules). Since the program was to be 
written in MS Pascal and FORTRAN* it could then be run on any 
IBM-class machine running the MS-DOS operating system* provided 
appropriate machine-specific graphic pixil and screen control 
module* were written. The success (and lack thereof) 
of this strategy is discussed below. 

During the fall of 1984# a Z-100 specific graphic pixel 
control routine was obtained from the Iowa State University 
Computation Center. Using this# Hogben wrote the module GRAPH, 
which sets up the graphics display and has procedures to graph 
points* lines* etc. Hogben wrote the module MENUS to display 
menus and interact with the user* and with advice from the 
Computation Center* wrote a brief Z^lOO screen control module 
As an illustration of designing a module and translating the 
module to programming language* Alexander wrote the module 
NULLCLN. In the testing of these modules* we discovered an error 
in the way double precision arithmetic was carried out in Pascal 
and Fortran in the version of MS--DOS we were using (MS-DOS v 1, 25 




Pascal and FORTRAN?? v. 3. 10). As there is no way to correct this 
errori and double precision is essential for the numerical 
methodSi it was necessary to obtain updated versions (MS-DOS 
V. 2. I3i Pascal and FORTRAN v. 3. 20). 

During the spring semester of 1985. Hogben and Tondra taught 
Math 51?X/ a graduate course covering numerical methods for 
studying phase planes and developing phase plane softuiare. 
Serious difficulties were encountered with both the quantity and 
quality of students. Only four students registered for the 
course (under a variety of credit arrangements) and two 
subsequently droppe; it. This created a serious shortage of 
manpower for writiig the modules, and impaired the design reviews 
(too few reviewers). 

Although th ? ablest student, James Coyle, dropped the class 
due to pressures of other commitments, he completed the main 
program before leaving, and participated in some of the later 
design reviews. There were difficulties with both the students 
who remained in the class. 

One of the students who remained in the class was receiving 
credit for a creative project required for his Master's degree, 
rather than course credit Although this arrangement worked well 
with a different student later* it led to delays in the project. 
When the student didn't complete his project (CRITPT) by the end 
of the course, he asked for more time (repeatedly). Although the 
underlying problem was that he lacked the ability to do the work# 
it wasn't until six months after the end of the course that he 
was finally relieved of his module and it was given to another 
student (Linda Ten Hoeve) as a creative Master'^i project. 
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The other student who remained in the course lacked 
essential programming skills. He uids given the job of adapting 
existing (public domain) numerical code for finding trajectories 
(Shampine 6..d Gordan* Computer Solution of Ordinary Differential 
Equations) for use in our software. Although he completed his 
uiorki it was so fTauied that it had to be essentially redone by 
Alexander later. 

Beacuse of the shortage of students* Hogben wrote the module 
for entering differential equations* PARSE. 

At the conclusion to Math 517x in M«»y 1985 the fol lowing 
modules had teen completed: MAIN, PARSE, GRAPH, MENUS. CRITPT 
was still assigned to the student who failed to complete it. 
During the summer of 1985* we hired Coyle to design and write the 
ARROWS module, which he did skillfully and efficiently. 

At the end of Fall semester 1985, the student assigned to 
CRITPT was removed from the project. Ten Hoeve was given the 
CRITPT module as a creative project required for her Master's 
degree. She completed this project successfully in the summer of 
1986. However, because of difficulties in working with the 
program in its entirity, she did not incorporate her procedures 
into the phase plane program, but wrote a small driver program. 

By the end of the summer of 1986, the program was complete 
except for the module PHASE and interfacing the Shamp ine/Gordan 
and Ten Hoeve subprograms into the main program. During the 
spring of 1987, Hogben adapted the Ten Hoeve code and Alexander 
adapted the Shamp ine/Gordan code. 

The only remaining part of the original design left undone 
was the PHASE module, to draw several trajectories. Since each 
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trajectory takes up to several minutes to draw, since the user 
can produce a phase plane of trajectories by repeating the 
trajectory command several timeSi and because the location of the 
trajectories desired depends on the differential equation, it uias 
decided to omit this feature. Thus the program itself uias 
complete (Z-lOO version). 

With the exception of the PHASE module, the final program 
follows the design developed during the fall of 1984. For 
technical reasons some modules are split between disk files and 
some disk files contain more than one module. 

Because of the difficulties encountered during Math517x. the 
students did not write the user manual as they went along^ as 
originally intended. Alexander is currently writing a detailed 
user manual for the program. When this is completed, we intend 
to test the software in an undergraduate differential equations 
course. In the meantime, brief instructions on how to use the 
program are available. 

In order to facilitate dissemination of the software, we 
intend to produce versions for other ns-^DOS machines, in 
particular the IBM-PC. The program design isolated machine 
specific features in two small modules, graphic pixel control and 
screen position control. During the summer of 1987. Hogben began 
work on the IBM-PC version of these modules. This was designed 
to be a minor job; but technical difficulties were encountered. 

The screen control module controls the position of the 
cursor on the screen, normal or reverse video, beeps, etc. All 
these features are simple to do on the Z-100. If the 
implementation of the IBM-PC had been similarly straightforward. 
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the IBM-PC version uiould have been quickly completed. Houieveri 
uie discovered that screen control on the IBM-PC is not 
automatically available. Although it can be arrangedi it does 
not uiork as well as the Z-lOOi and some features, e.g., reverse 
videoi are not* available on the graphics screen. We had been 
using reverse video to highlight the cursor position and error 
messages. The cursor position can be boxed with graphics and the 
error messages can appear normally, but separating out these two 
cases is requiring changes in other parts of the program. The 
IBM-PC version of the software should be completed this fall. 

The graphic pixel module for the Z-100 ujas obtained from the 
Iowa State University Computation Center. Fortunately, they also 
had available an IBM-PC version. This was assembled, tested, and 
is working well. 

Project Results 

The principal product of this project is the computer 
program, PHASE PLANE GRAPHICS FOR ORDINARY DIFFERENTIAL EQUATIONS. 
This program is currently available for the Zenith Z-100 and is 
being adapted for the IBM-PC and similar MS-DOS machines In 
order to obtain widespread distribution, uie are planning to 
distribute it commercially. It has been submitted to CONDUIT. 
This program is included in Appendiy B 

In the process of developing this program, we developed 
Denu format input/output and graphics control procedures callable 
from MS-DOS Pascal or FORTRAN These procedures may be useful to 
other software developers. They are included, along with a small 
demonstration program, in Appendix D 
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In the process of planning the coursei a variety of materials 

uiere developed. These include program design noteSi class plan, 

and a brief software development bibliography. These materials 

are included in Appendix C. 

Copies of the above materials may be obtained from 

Lesl ie Hogben 
Department of Mathematics 
louia State University 
Ames, lA 50011 

Tel. 515-294-8168 or 515-294-1752 
Because information about equipment is necessary to provide 
computer diskettes, such requests are more easily handled by 
telephon 

Summary and Conclusions 

This project successfully completed the development of 
planned softuiare, Phase Plane Graphics for Ordinary Differential 
Equations. Houiever, the intended method of development by a 
class of graduate students uias not completely successful. The 
class did not have enough students to write the program or 
provide effective design reviews. There were also technical 
difficulties encountered, with computer arithmetic and machine 
specificity. These difficulties were less serious. 

There a number of possible ways the student shortage could 
be overcome. Although information about the coursp was 
distributed to both the mathematics department and other 
departments, more personal recruiting through faculty advisors 
could have been done. Our best source of academic credit/unpaid 
student help was students who obtained credit for a Master's 
degree creative project, rather than actual course credit 
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There should have been more lead time before the course* to 
alloui It to be listed in the schedule of classes* and alloui 
students more time to plan to include it in thcur schedules. 

It is possible that uiith changes such as these that more 
students could have been obtained. Houieveri uie <=i>riously 
underestimated the rigidity of the graduate student program. 
Many students told us that they would like very much to take the 
coursei but they haj so many required courses that they couldn't 
fit it into their programs. 

In overcominj the student shortage difficulty, several 
promising methods for developing softuiare were utilized: 
Modifying currently available research progi^ams. having graduate 
students develop software as a creative project associated with a 
Master's degree, and hiring graduate students. 

Uell designed and thouroughly tested programs utilizing 
sophisticated numerical methods are widely available on mainframe 
computers, and are frequently in the public domain. 
Microcomputers have now become so powerful that many of these 
programs are accessible to them Although such programs 
frequently employ a high level of numerical sophistication, they 
range from difficult to impossible for an inexperienced person to 
use. Such programs are often in the form of subroutines that 
must be called from another program, thus rendering them 
inaccessible to the nonprogrammer. However, if combined with a 
user-^friendly program to ha^^dle input/output and a well written 
manua 1. existing numerical programs can provide the basis for 
valuable software. We successfully used the Shamp in e-Gor dan code 
for the finding of trajectories in Phase Plane Graphics for 
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Ordinary Differntial Equations. Another example of a research 
level numerical analysis program combined with a user-friendly 
input/output program is the MATLAD program for matrix operations 
using LINPACKi developed by Eugene Johnson and others at the 
University of Iowa. 

The ozhev strategy employed to successfully complete the 
software involved using graduate students to design and write 
parts of the program^ but not as part of the class. One student, 
Linda Ten Hoeve, wrcte the module for finding and classifying 
critical points as the creative project for her Master's degree. 
Another student, James Coyle, who had written the main program in 
the classi was hired to write the ARROW module. Both these parts 
were done well and on time, because the students had the 
necessary ability, programming skill, and sense of 
responsibility. Both these approaches provide good software 
development experience, albeit without the team approach and 
design review experience originally planned. Given the rigidity 
of graduate student schedules, it is unlikely that enough 
students could be recuited to work on Master's projects 
simultaneously to provide an effective group experience. It is 
possible that a group of students could be hired during the 
summer, although it is not clear that enough high caliber 
students with the requisite skills are available. This approach 
has been used by David Sharp and Colin Prowse, Southern Illinois 
University, in a FIPSE project involving computer environmental 
s imu lat i on. 

The results of this project suggest possible methods of 
developing undergraduate software by incorporating existing 
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research subroutines and using graduate students, either by 
hiring them or by awarding academic credit for individual 
projects. This project also developed a sophisticated and easy 
to use program for use in an undergraduate differential equations 
course. Phase PLane Graphics for Ordinary Differential Equations 
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student Produced Advanced Mathematical Softuiare 
Appendix A - Insights for FIPSE 

Although there uias substantial contact during the first year 
uiith FIPSE staff concerning FIPSE matters and the Technology 
Study Group/ thers was little contact concerning our specific 
project. There was also a lack of continuity of Program 
Officers, After Susan Forman lefti u^t? had contact uiith a variety 
of Program Officers. 

I'm not sure more contact about the project would have been 
particularly useful/ however. I do have a little concern that 
serious participation in the Technology Study Group could divert 
time and energy from projects. Although it is developing 
interesting ideas* our participation in FTSG was limited by lack 
of time available. 

The intrusion of other interests of the project coordinators 
was somewhat of a problem in this project All of us have a 
variety of interests and other projectsi which sometimes 
interfered with getting things done in timely fashion During 
the duration of this project, one of us wrote a book* one was 
assistant department chair, one was ov€»rloaded with graduate 
s>tudents, and two had leaves of absence 

One other remark about difficulties encountered. We found 
'•people problems'* such as lack of students much harder to deal 
with than technical problems. This may be due to partly to the 
backgrounds of the coordinators! who have more technical 
experience. 



student Produced Advanced Mathematical Software 

Appendix Q 

Phase Plane Graphics for Ordinary Differential Equations 

Contents 

Brief Instructions for Use 

Disk 1 

PHASEPLN EXE 
MACHINE. DAT 
HFLOOOOO. TXT 
HFLOOOOl. TXT 
HFL00002. TXT 
HFL00021. TXT 
HFL00003. TXT 
HFL00031. TXT 
HFL00032. TXT 
HFL00033. TXT 
HFL00034. TXT 
DEMO. EXE 

Disk 2 

MAIN. PAS 
PARSE. PAS 
□DENUM. PAS 
□DEFOR. FOR 
GRAPH. PAS 
MENUS. PAS 
SINGLE. FOR 
HSTART. FOR 
DIMACH FOR 
VNORM. FOR 
INTRP. FOR 
STEP2. FOR 
ZSCR. PAS 
ZGRAPH. ASM 
DEMO. PAS 
MENUS. DOC 
GRAPH DOC 



The source files ( PAS, . FOR, . ASM) can be read on any IBM-PC 
compatible microcomputer running MS-DOS The program itself can 
be run on a Zenith Z-100 running MS-DOS v 2 13 Place Disk 1 in 
the default drive and give the command 

PHASEPLN 



PHASE PLANt GRAPHICS FOR ORDINARY DIFFERENTIAL EQUATIONS 



L. Hogben; R Alexander, R Tondrai J Coyle; L. Ten Hoeve; and 

Iowa State Computation Center. 
This program was developed with support from the Fund for the 
Improvement of Post sec ondary Education, U S. Department of 
Education, and the Iowa State University Sciences and Humanities 
Research Institute All rights reserved, Iowa State University 
Research Foundation. 



SUMMARY OF INSTRUCTIONS FOR USE 

This program graphically displays solutions to the system of 
differential equations 

X ' ( t) = F( X, y ) y ' ( t ) G( x, y ) . 

INITIAL RUN 

If necessary, boot the system (MS-DOS v. 2 13). Insert the 
PHASE PLANE disk 1 in the default drive. Type 
PHASEPLN 

A title screen will appear, with the instruction ^at the bottom) 
to press RETURN. Do so, and thG graphics screen and main menu 
will appear. The top two-thirds of the screen is a region of the 
phase plane where things will be graphed The lower part of the 
screen is a menu that allows you to carry out various actions 
The last line supplies information, such as what you can do at 
this point or error messages 

You can move around the menu by using the arrow keys. To 
select an option, press the ENTER (or RETURN) key when the cursor 
IS on the desired menu item (as indicated by the fact that the 
this menu item appears in reverse video). You can obtain 
information about what various menu options mean by pressing the 
HELP key. 

To get the feel of the program, you should begin by moving 
the cursor to GRAPHICS MENU (press the down arrow key once), and 
select this option (press ENTER). This menu provides you with 5 
choices: ARROWS, SINGLE TRAJECTORY, NULL CLINCS, CRITICAL POINT 
and CLEAR GRAPHICS. Each of the first four provides another menu. 
Select SINGLE TRAJECTORY (by pressing the right arrow key once 
;:ind then pressing ENTER) This menu allows you to enter the 
starting point of the trajectory, X(0), Y(0), DRAW or ERASE the 
trajectory, enter the error parameters for the numerical 
routines, ABS ERROR, REL ERROR, or the duration of the 
trajectory, TIME T. 

De faults are supplied for all ciata The default 
differ en tial equation is 

x'(t) = (-'•x-y+4)/2 y'(t) = (x-y)/2 

Begin by selecting the DRAW TRAJECTORY option and watch the 
program draw the trajectory Note that while the program is 
graphing, the line at the bottom of the screen has changed It 
tells you it is graphing and ^^llows you to stop the graphing by 
pressing the ESC key (in case iz takes too long and you decide 
you don't want to wait for it to finish). 



Noui draui another trajectory by changing the initial point 
To do thisi you uiill have to select menu item X(0)* enter a neuj 
value (type in the number, ending with RETURN; see DATA ENTRY 
beloui for more information), select Y(0)> enter a neoi Vr»lut?. and 
iii?lGct DRAW again. If you want, you .nay also try to ERASE a 
trajectory* or change its length by changing TIME T 

Next* return to the graphics menu by pressing HOME This 
l<ey always returns you to the previous m^nu Then experiment 
with the other graphics options: ARROWS, k'hich allows you to draw 
the direction field of derivative vectors. NULL CLINE3, which 
allows you to uraw the null clines (x' O or y ' = 0), CRITICAL 
POINT* which allows you to find and classify a critical point 
(x' = O and y ' = 0) and mark it on the screen 

Next* clear the graphics screen* go back to the main menu* 
select SET DIFFERENTIAL EQUATION* and enter your own equations* 
by 5>electing X' = and typing in a new equation. (N. B 
multiplication must be represented by See below for a 

complete list). Do the same for Y', Now go back and try the 
various graphics options on your own equations. 

Finally* return to the main menu and select SET WINDOW. 
This provides a menu that allows you to modify the region of the 
phase plane that is graphed (as well as alter other features* 
^uch riS whether axes appear). Note that when you change the 
nrdr.jmum or maxixmum x or y values* the graphics screen does not 
change automatically. Once you have entered your new bounds* you 
must select REDRAW WINDOW to adjust the area graphed. This will 
al^o clear the graphics screen. 

When you have finished using the program* either turn off 
thf^ machine or return to the main program* select QUIT PROGRAM; 
and respond Y for yes when prompted 

P'urther information about the various menus is provided 
below* as is detailed information on entering data. 

DATA ENTRY 

You enter a piece of information (integer^ real number* 
(1 if f ewrential equation) by selecting a menu option that calls for 
the information and typing it in During data entry you cannot 
move around the menu or select other options. 

To enter a real number or integer* type in the number in the 
<jtu<:^1 form (not scientific n-^tation)/ ending ujith RETURN (or 
tINIER) If you make a mistake while typing <e g * a letter 
'nther than a digit)* the program will detect this and send you 
an error message to remove the error by pressing BASCKSPACE. If 
(before pressing RETURN) you decide that you do not want the new 
f^niry* you can press ESC and the program will revert to the its 
previous information (just prior to beginning entry). 

While entering a differential equation* you can move around 
in the equation by using the left or right arrows* delete the 
character left of the cursor with BACKSPACE (or DELETE)* delete 
the whole equation with DEL LINE* and insert characters by 
pressing the the appropriate keys. To end entry* press RETURN to 
accept your new equation or ESC to revert to the previous one A 
differential equation may include numbers (real or integer)* the 
v«iriables x and y* operands -»■* -***/,' * functions* par eritheses* 
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and the parameters P 1 1 P2, P3, P4; P5i P6 The table below lists 
symbols and their meanings The equation must also be in a form 
that the program can understand. It mill detect anything it 
cannot understanc as an error and supply an appr opr la'*"^ message 



bymb ol 
/ 

CXP{ ) 
LN( ) 
I 0Q{ ) 
SIN( ) 
C05( ) 
TAN( ) 
AT AN ( ) 
SQRT( ) 
ABS( ) 
( ) 

12 3 4 
6 7 8 9 



5 
0 



P1,P2,P3 
P4, P5i P6 



Mean ing 
add i t i on 

subtraction or negation 

mul t i p 1 i ca t i on 

d ivisi on 

exp onentiat i on 

exponential function 

natural logarithm 

base 10 1 ogari thm 

sine 

cosine 

tangent 

arc tangent 

square root 

absolute value 

parentheses (must be matched) 
digits for numbers 

decimal point 

blank (ignored except within number) 
var iab les 

parameters (values assigned by selectinq SET 
PARAMETERS) 



LIST OF MENUS 



MAIN MENU 

SET WINDOW. Menu for options controlling portion of phase 

p lane displayed 
SET DIFFERENTIAL EQUATION' menu for entering differential 

equa ti on 

GRAPHICS MENU, menu for graphics options 
CLEAR GRAPHICS, clears graphics screen 

QUIT PROGRAM quit programi requests c onf irrra t i on Y (yes) or 
N (no ) / def au It is no 

SET WINDOW Menu for options controlling portion of phase 
p lane d i sp lay ed 
MIN X = minimum x value graphed (left edge) 
MAX X =. maximum x value gi ^phed (right edge) 
X GRID =. number of grid points in the x direction 
MIN Y =• minimum y value graphed (bottom edge) 
MAX Y ~: maximum y value graphed (top edge) 
Y GRID = number of grid points m the y direction 
DRAW/ERASE AXES, draw axes if absent; erase axes if present 
REDRAW GRID: draui grid with X GRID and Y GRID points, 

does not clear grahics screen 
REDRAW WINDOW, draw window with indicated valuesi 
clears graphics screen 



SET 



DIFFERENTIAL E^iUATIQN 
SET PARAMETERS 



menu for entering differential 
equat 1 on 
menu for enter parmateters for 
differential equation 
enter equation for the derivative of x 
respect to time 

enter equation for the derivative of y 
respect to time 



uii th 



uii th 



SET PARAMETERS menu for entering parmateters for 







d i f f erent lal 


equat i on 


PI =: 


enter 


parameter 


PI 


(real 


numb er ) 


P2 


enter 


parameter 


P2 


(real 


numb er ) 


P3 =: 


enter 


parameter 


P3 


(real 


numb er ) 


P4 - 


enter 


parameter 


P4 


(real 


numb er ) 


P5 =: 


enter 


parameter 


P5 


(real 


numb er ) 


P6 =: 


enter 


parameter 


P6 


(real 


numb er ) 



GRAPHICS MENU, menu for graphics options 

ARROWS: menu for drawing direction field of derivative vectors 
SINGLE TRAJECTORY: menu for drawing single trajectory 
NULL CLINES: menu for drawing null clines 

CRITICAL POINT; menu for finding and drawing critical point 
CLEAR GRAPHICS: clears graphics screen 



ARROWS, menu for drawing direction field of derivative vectors 
ARROW DENSITY enter integer to control arrow density; 

larger means more smaller arrows 
DRAW ARROWS: finds and draws direction field 
ERASE ARROWS, finds and erases direction field 



SINGLE TRAJECTORY: menu for drawing single trajectory 

X(0) =: enter the X "~c oor d 1 na t e of the initial point (real number) 
Y(0) =. enter the y-coordinate of the initial point (real number) 
DRAW TRAJECTORY find draw path from t=0 to t=:TIME. 
ABS ERROR = error tolerance for trajectory finder (real 

number; at least one error tolerance must be >0). 
REL ERROR =: error tolerance for trajectory finder (real 

number! at least one error tolerance must be >0). 
ERASE TRAJECTORY, find erase path from t=0 to t=^TIME 
TIME T =• end time for trajectory 



NULL. CLINES: menu for drawing null clines 

NUMERICAL SEARCH SPACING enter space between points tested for 
sign change in the derivative (integer* larger - more 
space between marks* fewer* less accurate marks) 

DRAW NULL CLINES. finds and marks null clines. 

ERASE NULL CLINES. finds and erases null clines. 
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CRITICAL POINT, menu for finding and drau/ing critical point 

NEARBY X =. enter x value of point near critica? point sought, 

'•NEARBY" changes to "FOUND'' and value changes after 

the critical point is found 
NEARBY Y «=: enter y value of point near critical point sought/ 

"NEARBY" changes to "FOUND" and value changes after 

the critical point is found 
DRAW CRITICAL PT: draw critical point (must be found first) 
ERASE CRITICAL PT. erase critical point (must be found first) 
FIND CRITICAL PT. find and classify critical point; changes 

"NEARBY'' to "FOUND". 
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student Produced Advanced Mr>themat ical Software 

Appendix C 
Material? Developed for Classroom Use 

Contents 

Software Design Notes 
Software Review Checklist 
Class Schedule 

Software Design Bibliography 
Design Notes for Module NULLCLN 
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General Requirements for Math 517X (FIPSE Project) Software 

L. Hogb en January 1985 

This document was developed with support from the Fund for the 
Improvement of Pos tsec ondary Education. U. S. Department of 
Educationi and the Iowa State University Sciences and Humanities 
Research Institute. All rights reserved/ Iowa State University 
Research Foundation. 



Spec if ications 

Specifications will be provided by the 
detailing the modular structure of the project. For 
those procedures intended to be available to other 
programs will have purposeSi declarations^ calling 
parameters specified 



inst r uc t or s^ 
each module> 
modules or 
syntaxi and 



Oes i gns 

Design will be done during the course. Designs will b 
adequately documented; following the style of the procedure 
NULCLN. Prior to coding the design^ a design review will be 
conducted by the class After the design is accepted it will be 
coded. 



Code 

All code will be written in MS-PASCAL or MS-FORTRAN and will 
be capable of being run successfully on any MS-DOS machine with 
192K or more memory. Should it be necessary to provide machine 
specific information* this will be isolated in a datafile <e. g 
machine constants) or a minimal machine-specific code module 
(e.g. ZSCREEN); commonly called a device handler. 

All code will conform to the layout and comment style 
exhibited in modules MENUS and GRAPH. Specifically, most 

CO mm en ts will be one line exp lanat ions of a particular action 
interspersed in the code* except that an initial explanation of 
all global module variables will be given. Each level of 
dependence will be indicated by indenting 2 spaces. For example, 
under an IF-THEN statement, indent an additional 2 spaces. 
Separate statemnts will be on separate lines. 

After the code is written, grammar will be tested by 
comp iling it and any grammatical flaws will be corrected. 
Compilable code will then be presented presented to the class for 
a code review. During the code review, structural and functional 
testing will be designed. 

Test ing 

After the code is accepted it will be tested. Any errors 
found during testing will be corrected and the module retested 
On completion, the software will be tested as a whole. 



User Interface 

Material will be presented to the user by a series of menus. 
Each menu will appear at the bottom of the screen with a command 
line. The upper portion of the screen will be a graphics area. 

The software must be crashproof. In case of error it should 
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provide helpful diagnostic messages Program accessible 

helpfiles uiill be be part of the softuiare package. AnM action 
the user takes should be reversible. An exception to this 
requirement (e.g. terminating the program) should reuire user 
confirmation. All I/O must be done in such a manner as to 
prevent screen scrolling. 

A clear« we 1 1-uir i t ten user manual uiill be written. It will 
explain to the user the function of each menu and commandi 
illustrating these functions with examples. 

Modular Structure 

The main prop'^im will be responsible for initializing the 
menus and machine constants. It will also be responsible for 
controlling the transition between menus. The authors of the 
ma in program will be responsible for providing an appropriate 
helpfile for each menu. The ma in program must be written in 
Pascal. It will have access to modules to cirry out the 
graphics^ menu layout/ and necessary ma t hema t ic s. 

Each module will be responsible for its own error handling 
and diagnostic messages. 
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REVIEW CHECKLIST 



Note anything you do not understand. 

Note anything that is not sufficiently well defined to enable you 
to carry out the next stage or to check with *"he previous stage 

Does it conform to tht general standards? 

Does it conform to previous stagesT-' 

Are theT»e any unstated assumptions or restrictions? 

Is it well structured? 

Is it testable? 

Is it maintainable? 

I/O: 

Is it crash proof? 

Are useful messages supplied to the user to assist in 

recovery from errors? 
Is it flexible? 

A^i thmet ic Errors: 

Is it crashproof? 

N'Jmer ica 1 Algorithms: 

Is it crashproof (i.e guaranteed correct to within some 
tolerance)? 

Are useful messages supplied to the user to assist in 
recovery from errors? 

Is It hardware compatible'^ 

Is it software compatible? 
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MATHEMATICS 517X 



Ho^beivTandia 



QA76.6 
J66 

QA297 

M527 

1984 

QA76.6 

D845 

1984 

QA76.6 
F86 

QA372 
S416 



Jones - Software development, a rigorous approach 



Miller - The engineering of numerical software 



Dunn - Software defect removal 



Wulf - Fundamental structures of computer science 



Shampine - Computer solution of ordinary differential equations 
The initial value problem 
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Dnft doig? for NULCLN procedure to ploc nuUcUiiei 



Dmffi Review ComiDenCi 



ILK. Alennder 



January 1965 



(Hogben) 



Thii document wm developed with support from the Fund for the 
Improvement of Pottaeooodaiy Education, U. S. Department of 
Education, and the kwa Sute UnivcrBity Sdencet and Humanitiea 
Rcaearch Institute. All righu rcMrved, Iowa Sute Univeiaity 
Research Foundation. 



on all 

real *> im conv 
(including int/mt) specify 
rounding 



OVERALL STRATEGY 

Search for sign changes between the pouiu of a square grid placed on 
the window. When a sign change is found in either dertvathre, localize it 
fiirther by a single secant step and mark it by three pixels in a vertical 
column if diAlt « 0 or three pixels in a horizontal row if dy/dt » 0. 



IMPORTS 



XO,YO,Xl,Yl: REAL; 
SPACE: INTEGER; 



{Cborinates of window pCO,Xl] x [YO,Yl] } 

{Size in pixeb of grid on which changes } 

{of sign are sought; user selecu one of } 

{4,6,12,16,20 from menu. } 



what does secant 

step mean - 

cf. seekz lin interp. 



does not 
conform to 
main specs 



EXTERNAL and FUNCTIONAL REFERENCES: 

DERIV {Subroutine to evaluate derivaives} 
GRNON {Procedure to turn on uxUvidual pbteU} 
SGN {Signum function} 

MIN {INTEGER minimum of INTEGER argumenU} 
LOCAL IDENTIFIERS: 



VECTOR means ARRAY[1..2] OF REALS; 

X: VECTOR; {Position of current point in grid} 
DX: VECTOR; {Vector of distances in the coordinate directions } 
{ } 
{DXIl] - ((Xl-XOyMO) • SPACE } 
{DXI21 - ((YJ-YOyiSO) • SPACE } 
XDOT: VECTOR; {Derivatives evalutcd at current point } 
OLDRV: VECTOR {Derivatives evaluated at 'north' neighbor } 

{of current point } 
SAVDRV: ARRAY[0..45] OF VECTOR {Derivathws at grid poinU in the 

{next oolunm to the left It is 
{assumed tha he window will occupy at most pixels 0 to 
{179 in the vtidcal direction so that all derivatives in a 
{column of grid points with the densest spacing (4) can be 
{saved in this array. 
RTEDGE: INTEGER CONSTANT, {Rightmost pixel column « 639} 
BMEDGE: INTEGER CONSTANT, {Bottom pixel row - 179} 
IX,IY: INTEGER; {Counters for pixel position in directions } 
INDY: INTEGER; {The ordinal number of the current grid point in } 
{its column, counting the top point as C. Used } 
{as an index into the amy of saved derivatives } 



Machine constants 
should be passed 
this should be in body 



Good note of restriction 
other machines 



Import this 
Specify integer type 
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c<0 



LOCAL PROCEDURE: 



SEEKZ(XJ>X,DUYJ>IR,XDOT,OLDRV); 

{ Checks each oompoiient of the derivative for a sign change in the } 

{ directioo DIR (-1 or 2). If one it found, localize it further by } 

{ a lecant step, and turn oo the appropriate pixels on the soeen. } 

{ Details for SEEKZ are ghm after design of NUCLN. } 

BEGIN { NULCLN } 

{ Initialize DX } of previous page 

{ } 

{ Loop through the columns in the grid ) 

IX:-Q; X(1]:«X0; use separate lines 

WHILE DC < RIBDOE DO { Loop through the grid poinu down this column } use For instead 

lY 0; X12] Yl; INDY 0; 
WHILE lY < BMEDGE DO 

DERIV(X^AR,XDOi:); { Evaluate derivathws at new grid point 
{ 

IF lY > 0 THEN {There is a grid pt to the north; check 
SEEKZ(X4>X,IX,IYAXDOT, OLDRV); { for sign changes. 

{ 

{ Check for sign changes to the left of the current point 
{ 

IF DC> 0 THEN SEEKZ(XJ>X,DCIY,l,XDOT, SAVDRV{INDY]); 

( } 
{ Save the derivatives bom this point: } 
OLDRV <- XDOTi 
SAVDRV[INDYJ <~ XDOT, 
{ > 
{ On to the next point in the column } 
lY lY + SPACE; 
INDY :- INDY + 1; 
XI2] :- X(2] - DX(2]; 
END WHILE { lY < BMEDGE } 
{ } 
{ On to the next column } 
DC :« DC + SPACE; 
XII] :- XII] + DXIl]; 
END WHILE { DC < RTEDGE > 
END NULCLN 



INTERNAL PROCEDURE SEEKZ: 

X4>X,XDOT,OLDRV,DCIY as in NULCLN 

DIR: INTEGER; { Index of coordinate search direction » 1 or 2 } 



LOCAL IDENTIFIERS : 

K: INTEGER; {Index on the oomponenU of the derivative } 
PIXEL ARRAY(L^] OF INTEGER: { Coordinates of pixel to light on screen } 
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BEGIN (SEEKZ) 

FOR K :- 1 TO 2 DO { Check tuh component } 

IF SGN(XDOT[iq) < > SGN(OLDRV[K]) THEN { Sign change found } 



PIXEL(DIRJ PIXEL(DIR] - OlV^^'^> XDOTfKI # SPACE 

XDOTpq -I- OLDRVpq 

{ The power of (4) it became Y decrentes at you go down the tcreen 
{ This fonnula amounts to a secant step 

{ Hie pixd in the other direction should be t'je common DC or lY 
( index of the grid pointt being checks. 

{ . 

{ If X " 0, illumine 3 pixels in a vertical column centered here. 
{ . 

{ If y - 0, turn on 3 piieli in horizontal rown centered here. 
{ "Here* means at location (P1XEL{1]4'IXEL(2]) on the 
{ 



incorrect 

unclear; linear interp 
of locof 0 
TRUNC ( +0.! 



SumiDuy of chiniei to Gnt dnfl detign oC NULCLN plotter (Alannder) 



Importi: Hie pixel boundariet of the window hive beeo added. I don't 
MMune any more that the window to [0,639] x [0,179]. 

The dnw/enae flag ia ^dded. The fonner detign drew only. 

The vector PAR of mer-aet panmcten in the derivative functiona 
ia needed, of oouiae. 

Body; DX(1], DXf 21 are now defined to be the acndni of tingle wxela 



SEEKZ pnxx Delete tpecutor paramelen X, DX 

Add grid tpadng, pixel #8 of window boundarict, draw^ite flag. 



In addition to recognizing lign changet, the detign now handlea the 
caae of a derivative being aero at both grid pointa. The detign 
it rather crude - thit rare cate it tetted finit; and interior grid 
poinu at which a derivative ia exactly aero wiQ be mariced twice. 

The formula for the lecant ttep wat incorrect; it it now corrct 

MARK proc Thk new, lower-level routine hat been added to actually perform 
the drwAerate on the tcreen. !t utet MIN A MAX to enture that 
GRNON A GRNOFF are only called with pixel addrettet inside the 
window. Detailt are given which did not appear in the fittt drha. 



along the coordinate aiea. 
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Design Revi<^; Commentt 



(Hogbcn) 

Memo: Mathenuitks 517X 

From: Roger Akannder 

Re: Revised Design of NuUdine Procedure 



OVERALL DESCRIPTION: The user has selected a grid sfwctng of 

4A12,16 or 20 pixds finom Che meou. This routine searches - why only 

for sign changes in both derivatives between points of a grid 4, 8, etc 

with the selected spacing. When a sign change is found, the 

root is further localized by one secant step, and then marked by 

three piids in a vertical column if dx^t - 0, or three pixels 

hi horizootal row if dyAh » 0. The pixds are turned on if 

the user has requested "DRAW, turned off if the user requested 

"ERASE". 



IMPORTS: 

X0,X1,Y0,Y1: REAL; { Coordinates of window pCOPCl] x ) 

{ |Y0,Y1) ) 
LFEDGEJITEDGE3MEDGE,TPEDGE: INTEGER: 

{ Pixel numben corrcspondhig to X0,X1 

{ and Y0,Y1 respectivdy. 
SPACE: INTEGER { Grid spacing m 4^12,16 or 20. 
DRAW: BOOLEAN { True -> Draw; False -> Erase. 
PAR: VECrOR(6) of REALS; {Parameters hi derivatives. 

EXTERNAL AND FUNCTION REFERENCES: 
EVAL {Subroutine to evaluate dertvait^cs ) 
GRNON,GRNOFF { Pnxsedures to turn off {resp. on) single pixels 
SGN { Signum function 

MIN { INTEGER minimum of INTEGER argumenU 

DBLE { INTEGER or REAL to DOUBLE PRECISION conver ion 

LOCAL IDENTIFIERS: 

VECTOR means ARRAY[L.2] of DOUBLE PRECISION 



good 

use 4 < SPACE < 20 



turn on (resp. off) 



X- VECTOR; 
DX: VECTOR; 

XDOT: VECTOR; 
OLDRV: VECTOR; 



{ Position of current grid point 
{ Vector of disunoes hi the coordinate directions 
{ represented fay one pixel 
{ Derivativts evaluated at current point 
{ Derivatives at 'north' ndghbor of X 
SAVDRV: ARRAY[0..45] of VECTOR; 

{ Saved derivath«s at grid pomts hi the pre^ous 
{ oolunm to the left The size of this array is 
{ based on the assumption that the window 
{ occupies at most pittls 0 to 179 m the veitical 
{ direction. 

{ Pixd position in directions. 
{ The oidbal number hi iu column of the current 
{ grid pohit, counting the top pouit as 0. 
{ Used as an fakkx hito the array of saved 
{ derivatives. 
MINPX, MAXPX: ARRAYIL.2) of INTEGER; 

{ Pixel boundaries of window ) 



DC,IY: INTEGER; 
INDY: INTEGER; 
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begin { NULCLN ) 






{ Conpule X- and y- ipndng of padt) 






DXll] (XI - XO) / (RTCDGE - LFEDOE); 






DX[21 (Yl - YO) / ^MEDOE - T7EDGE); 






MINPX[11 :- LFEDGE; 






MINPX(21 T7EDOE; 






MAXPXPI RTEDGE; 






MAXPXPI BMEDOE; 






{ Loop through oolumm of the grid ) 






for DC LFEDOE to R1EOOE hi steps of SPACE do 






{ Compute x-coonUnate of grid poinU in this column ) 






X[ll XO + DX[11 • (DC - LFEDOE) ; 






{ Loop through grid pohits in this column ) 






for lY T7EDOE to BMEDOE fai steps of SPACE do 






{ Compute hidcat and y-axxdfaiate of this pobt ) 






INDY (lY - T7EDOE) / SPACE; 






X[21 Yl - DBLE(IY - T7ED0E) • DXp]; 






EVAL(XJ>AR,XD01> {Bvahuite derivath^ at this pt ) 






if lY > T7EDOE then { Check for sign-changes to north ) 






SEEKZ(DC,IY^PACEJ4INPX,MAXPX,2,XDOT,OLDRV,DRAW); 






fi; 


fi^endif 




if DC> LFEDOE then { Check for sign<hanges to the left ) 






SEEKZ(DCIY^PACE>fINPXJ4AXPX,l,XDOT, 






SAVDRV(INDY]4>RAW); 






fi; 






{ Save derivatives from this pohit ) 






OLDRV <- XDOP, SAVDRV(INDY] <~ XDOP, 






od;{IY) 






od;{DC) 






end; {NULCLN ) 






proc SEEKZ (DC,IY: INTEOER; { Pixel cooidinates of current pt ) 






SPACE: INTEGER; { Grid spacing ) 






MINPX, MAXPX: ARRAY[L^1 of INIBOER; 






{ Pixel boundaries of window ) 






DIR: INTEGER; { Coofdinate direction in which ) 






{ to search: 1 or 2 ) 






XDOT.OLDRV: VECTOR { Derivatives at current grid ) 






{ point and previous point } 






DRAW: BOOLEAN); { Draw/ersse request ) 






local identifien: 






K: INTEGER; { Index for componenU of derivative ) 






PDCEL: ARRAYIL.2] of INTEGER; 






{ Screen coordinates of zero of a derivative ) 
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ERLC 







beiin { SEEKZ) 

forK:-llo2do { Check each component ot derivative } 



PKELfl] :- DC; PIXEL(2] :- lY; 

if XDOT[K]«0 and OLDRV(K]-0 Uien { Tfto zeros ) 

MARIC(PIXEUMINPXMAXPXJU>RAW); 

PDCELpiR] PIXEL[DIR] • SPACE; 

{ Back up PSACE ptadt in the DIR direction ) 

MAR]C(PDCEUMINPX>fAXPX^RAW); why do it twice? 

ebe 

if SGN(XDOT[iq <> SGN(OLDRV(iq then 

{ Sip chanfe found in thii component Back up ) 
{ by a secant step parallel to the DIR axis. ) 
PDCELfDlR] :- PIXEL(DIR] - SPACE • 

XDOTpq / (XDOTpq • OLDRV[Kl); 
MAR]C(PIXEL»MINPXJ^fAXPX>fCJ)RAW); 
fi; { Sign change ) 
G; { zeros } 
od; { Check components ) 
end; { SEEKZ) 



proc MAR1C(PIXEL: ARRAY(1..2] of INTEGER; 

{ Screen location of zero of derivative } 
MINPXAfAXPX: ARRAY(1.^] of INTCGER; 

{ Boundaries of wmdow on screen ) 
K: INTEGER; { Index of vanishing derivathv component ) 
DRAW: BOOLEAN); { Diaw/erase request ) 



local idcntifierK 

PFDC ARRAY[U2] of INTEGER; { Pixels to alter ) 

L: INTEGER; { Loop index to alter 3 pixels ) 

NOTX: INTEGER; { The index (either 1 or 2) which is not K. ) 

{ If XDOTpq is zero, the vectorfield ) 
{ is perpendicuhu- to the K coordinate axis. ) 



begin (MARK) 

PFTXpq :- PDCELpCJ; 
NOTX :- 3 - K; 

for L -1 to 1 do { Fk 3 pix. centered at the zero ) 

PFTXJNOTK] :- PIXEL(NO + L; 

{ But don't get out of the window ) 

PFDCJNOTK] :- MIN(MAXPXINOTK]J»FTX(NOTXl); 

PFTXJNOTK] :- MAX(MINPX|NOTK]J»FDC[NOTK]); 

if DRAW then GRNON(PnX(lJ,PFIX{2J); 
else GRNOFF(PFIXIll,PFlX(2J); 

il; { DRAW ) 
od; (L) 
end; { MARK ) 



ERLC 
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Remirki on 13 JAN 85 dnft design of nuUdine plotter. (Aknnder) 



I. There it no cheddng of aifumentt. Pottibilitiet include: 

a) create a type, SPACETYP with only values 4^12,16, A 20. 
Then the compiler can check SPACE. 

b) type XPDCEL - (0..639]; YPDCEL - [0.224]; 

hawe the routine actually check LFEDGE < RTEDGE, 

TPEDGE > DMEDGE 



2. PASCALrlike decUrationt of edemal A function referencei have not 
been given. It ia likely that this ia desirable for checking argument 
types. Some possible pix)blenis with functioanlity, toa* 

a) SGN muat implement the mathematical, 3-vdued signum function. 
Fortran's SIGN (which takes two argumenta) may not do the trick. 

b) The design is completely stoppy about REAL4/REAL8 and INTEGER2/INTEGER4. 
What INTEGER" meana depends on the implemenUtion language! the MACRO- 
assembler code for GRNON gives the calling sequence as 



I would be veiy surpriied if both worked! Conclusion: we need a 
language-independent way to specify argument types in "absolute^ terms, 
a.;d we must pay careful attention to type conversion at the design 
stage, or the code stage will be a killer. 



XO < XI 
YO < Yl 



FORTKAN: 
PASCAL: 



CaU GRNON(IX,IY) 
GRNON(X,Y: INTEGER4); 



This is wrong 



ERIC 
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